Numerical Simulation for Porous Medium Equation by Local Discontinuous Galerkin Finite Element Method

In this paper we will consider the simulation of the local discontinuous Galerkin (LDG) finite element method for the porous medium equation (PME), where we present an additional nonnegativity preserving limiter to satisfy the physical nature of the PME. We also prove for the discontinuous P0 finite element that the average in each cell of the LDG solution for the PME maintains nonnegativity if the initial solution is nonnegative within some restriction for the flux’s parameter. Finally, numerical results are given to show the advantage of the LDG method for the simulation of the PME, in its capability to capture accurately sharp interfaces without oscillation.